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Recent data on charmonium production in B-meson decays suggest that charmonium hybrid 
mesons with mass ~4 GeV may be produced in B-decay via cc colour octet operators. Some 
of these states are likely to be narrow with clean signatures to J/ipir + Tv~ and J/tpJ final states. 
Experimental signatures and search strategies for existing B-factories are described. 
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I. INTRODUCTION 



The existence of gluonic excitations in the hadron spec- 
trum is one of the most important unanswered questions 
in hadron physics |l| . Hybrid mesons form one such class 
which consists of a qq with an excited gluonic degree of 
freedom. Although there is mounting evidence for hy- 
brids consisting of light quarks they still await confirma- 
tion 0, H, HI ■ Recent observations of charmonium states 
in exclusive .B-meson decays H S S S S E3 indicate 
that charmonium is produced significantly in the color 8 
leading us to argue that charmonium hybrids (i/y ) pdj 
should also be produced in B- meson decay The 
unambiguous discovery of such a state would herald an 
important breakthrough in hadronic physics, and indeed, 
in our understanding of Quantum Chromodynamics, the 
theory of the strong interactions, ft would also provide 
important input to refine models of hadron structure. 
In this letter we examine the production of charmonium 
hybrids and how one might observe these new types of 
mesons in i?-decays. While various elements of our ar- 
guments have appeared elsewhere, by putting all the in- 
gredients together, a more complete picture of charmo- 
nium hybrid production and detection emerges that, we 
hope, will encourage experimenters to pursue the neces- 
sary analysis. In some cases we can only resort to what 
are at best order of magnitude estimates of various pro- 
cesses. However, we expect that they will provide useful 
guidance for the initial exploration of this new physics 
frontier. Our goal is to point out that it may be possi- 
ble to discover charmonium hybrids in B decay and to 
suggest likely signatures to do so. We review their spec- 



troscopy, argue for their production in B-decays, and sug- 
gest experimental strategies for detecting charmonium 
hybrid mesons. 



II. HYBRID CHARMONIUM SPECTROSCOPY 

Lattice gauge theory and hadron models predict a rich 
spectroscopy of charmonium hybrid mesons [ill IhH ITU 
[3 113 Ell! Ill EI- For example, the flux tube model 
predicts 8 low lying hybrid states in the 4 to 4.2 GeV 
mass region with J PC = ±=F , 1 ±=F , 2 ±= F, and 1 ±± . Of 
these states the O" 1 , 1 and 2 H have exotic quan- 
tum numbers; quantum numbers not consistent with the 
constituent quark model. The flux-tube model predicts 
M{$ g ) -4-4.2 GeV 0IH; lattice Q CD predictions 
for the J PC = 1 -+ state range from 4.04 GeV to 4.4 GeV 
01 with a recent quenched lattice QCD calculation 
H3 finding M{\-+) = 4.428 ± 0.041 GeV. These results 
have the 1 h lying in the vicinity of the D**D thresh- 
old of 4.287 GeV. There is the tantalising possibility that 
the I h could lie below D**D threshold and therefore be 
relatively narrow. 



III. HYBRID PRODUCTION 

Recent developments in both theory and experiment 
lead us to expect that charmonium hybrids will be pro- 
duced in B decays. The partial widths for B — > cc + X, 
with cc representing specific final states such as J/ip, ip', 
Xco, Xcij Xc2, 3 £>2, 1 -£?2 etc., have been calculated in the 
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NRQCD formalism El El El El El E3 which factor- 
izes the decay mechanism into short (hard) and nonper- 
turbative (soft) contributions. The hard contributions 
are fairly well understood but the soft contributions, in- 
cluded as colour singlet and colour octet matrix elements, 
have model dependent uncertainties. Insofar as hybrid 
cc wavefunctions have a non-trivial colour representation 
they can be produced via a colour octet intermediate 
state. 

Over the last decade there has been great theoretical 
progress in the understanding of exclusive B decay to hid- 
den charm final states El El El El El El E3 Central 
to this is the recognition of the importance of the colour 
octet contributions to these decays. Although the colour 
octet terms are higher order in the velocity expansion 
for the soft contributions, the Wilson coefficients for the 
colour octet subprocess are significantly larger than that 
of the colour singlet subprocess in the hard contributions 
to the decay. The net result is that the colour octet com- 
ponents play an important role in these decays. This has 
been dramatically confirmed by the observation by the 
CLEO, Babar, and Belle collaborations [1I1I1HHE3 
of the decays b — » Xco and b — > \c2 which proceed via 
colour octet and have been measured to have compara- 
ble branching fractions to the decays b — > ip + X and 
b — ► Xci which have sizeable colour singlet components. 
Calculations of the BR's to the \D(cc) states, which are 
higher order in the NRQCD expansion but have a large 
colour octet contribution, also predict BR's roughly com- 
parable to those of other charmonium states [Hi El E3 • 

The NRQCD approach assumes a Fock space expan- 
sion of the state which gives a velocity scaling of the 
various wavefunction terms contributing to the soft pro- 
cess. For conventional mesons the colour singlet term is 
the leading term in the expansion and the colour octet 
term is higher order in d. In contrast, in many models 
describing hybrid states with constituent glue, such as 
the bag model [3 and constituent gluon models [2lj |. 
the colour octet cc configuration is the leading term. In 
the flux tube model the factorization of the QQ colour 
is not so clear and the main uncertainty is estimating 
the hadronic matrix elements. Chiladze et al. [l3| esti- 
mated the octet matrix element by explicitely calculat- 
ing the wavefunction in an approximation to the Isgur- 
Paton flux tube model 01 while Close E3 estimated 
the S-wave matrix element by rescaling the P-wave octet 
matrix element by 1/v. Chiladze et al. |l3j estimated 
the branching ratio B\B -> ip g (0 + ~) + X] ~ 1CT 3 for 
M ~ 4 GeV (though recent quenched lattice calcula- 
tions suggest M(0+-) = 4.70 ± 0.17 GeV, and hence 
will be inaccessible). Note that the decay to the H 
hybrid is suppressed by a spin factor of g while hybrid 
states with higher spin have larger statistical factors lead- 
ing to larger branching ratios. Close et al. ^2 es ti- 
mate a similar branching ratio to 1 h and argued that if 
M g < 4.7 GeV, the total branching ratio to tp g for all J p 
could be B[ip g (aH J PC ) + X] ~ 0(1%). Thus, using two 
different approaches for estimating B[B — > ip g + X] both 



Chiladze et al. 13] and Close et al. [l2| obtain similar 
results. Both calculations estimate BR's of 0(0.1 — 1%) 
which are comparable to the BR's for conventional cc 
states. We conclude that charmonium hybrids should be 
expected to be produced with roughly the same branch- 
ing fractions as conventional charmonium states. 



IV. DECAYS 

There are three important decay modes for charmo- 
nium hybrids: (i) the Z weig allowed fall-apart mode 
xp g -> £>(*>**)!)(*•**) [H E2 E2; (ii) the cascade to 
conventional cc states, of the type ip g — ► (cc)(gg) — > 
(cc) + (light hadrons) and i/j g — > (cc) + 7 E3 > ( m ) decays 
to light hadrons via intermediate gluons, ip g — > (ng) — > 
light hadrons, analogous to J/ip — > light hadrons and 
rj c — > light hadrons. Each mode plays a unique role. ip g 
hybrids with exotic J PC quantum numbers offer the most 
unambiguous signal since they do not mix with conven- 
tional quarkonia. 

(i) Decays to D^D^*h In addition to J PC selection 
rules (for example, 2 h and 2~~ decay to DD are for- 
bidden by parity and the exotic hybrid ip g (0 ) decays 
to final states are forbidden by P and/or C 
conservation) a general feature of most models of hybrid 
meson decay is that decays to two mesons with the same 
spatial wave function are suppressed 33]. Therefore, de- 
cays to DD should not occur, but small differences in 
wavefunctions could lead to small but finite widths to 
DD*. Seeing DD* but not DD or D*D* would be a 
striking signal for a hybrid meson. The dominant cou- 
pling of charmonium hybrids is to excited states, in par- 
ticular £)(*)(£ = 0) + D**(L = 1) states for which the 
threshold is ~ 4.3 GeV. This is at the kinematic limit 
for most mass predictions so that decays into the pre- 
ferred £)(*)£)** states are expected to be significantly 
suppressed if not outright kinematically forbidden. Es- 
timates for the various ipg decay widths and branching 
ratios are given in Table I . In columns 1-4 a hybrid 
mass of 4.1 GeV is assumed, which is below DD** thresh- 
old, while column 5 and 6 employs a mass of 4.4 GeV, 
above DD** threshold. This enables us to gauge the 
model dependence of the results and the effect of open- 
ing up the DD** channel on the total width. The original 
Isgur Kokoski Paton flux tube model [3(j predicts partial 
widths of ~ 1 — 20 MeV, depending on the J PC of the 
hybrid |3lj] while a refined version of this model predicts 
smaller partial widths of 0.3-1.5 MeV These widths 
are quite narrow for charmonia of such high mass. If the 
hybrid masses are above D** threshold then the total 
widths increase to 4-40 MeV for 4.4 GeV charmonium 
hybrids which are still relatively narrow. The challenge 
is to identify decay modes that can be reconstructed by 
experiment. 

(ii) Decays to (cc) + (light hadrons): The ip g — > (cc) + 
(light hadrons) mode offers the cleanest signature for ip g 
observation if its branching ratio is large enough. In ad- 
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dition, a small total width also offers the possibility that 
the radiative branching ratios into J/tp, rj c , XcJi an d h c 
could be significant and offer a clean signal for the detec- 
tion of these states. 

For masses below DD** threshold the cascade de- 
cays ip g — ► (ip, T) c , . . .) + (gg) and annihilation decays 
tp g (C = +) — > (gg) —> light hadrons will dominate. If 
the masses of exotic J PC states are above DD** thresh- 
old their widths are also expected to be relatively nar- 
row for states of such high mass, in which case cas- 
cades to conventional cc states transitions of the type 
ip g —¥ (ip, ip') + (light hadrons) should have significant 
branching ratios |29j making them important signals to 
look for in tp g searches. The hadronic transition rates to 
conventional charmonia from either C = + or C = — 
are similar because both are the same order in a s and 
the charge conjugation of the conventional cc daughter 
should be the same as that of the cc hybrid parent since 
two gluons (C = +) arc emitted in the lowest order pro- 
cess. 

A transition between two quarkonium states proceeds 
via the emission of gluons by the heavy quark and the 
subsequent conversion of the gluons into light hadrons 
[34|. The emission of the gluons is typically treated as 
a multipole expansion of the colour gauge field to esti- 
mate rates for hadronic transitions between QQ states. 
Kuang and Yan j3^ estimated the matrix elements be- 
tween quarkonium states by inserting intermediate states 
with the string in its vibrationally excited lowest mode, 
ie. hybrid states. Thus, the matrix elements for hadronic 
transitions between conventional quarkonia are related 
to hybrid-conventional quarkonium hadronic transitions. 
The widths (cc) -> nnJ/ip are typically O(10- 100) keV 
while (cc) -> rjJ/^P are typically O(10) keV H. (BES 
has recently measured r(V>(3770) -> J/ipn+TT-) = (139± 
61 ± 41) keV Hi] which is a 3 L>i S^ir transi- 

tion.) It seems reasonable to assume that the partial 
widths for the decays ip g (l — > r/ c + (irir, r), rf) and 
ip g (0^ , 2" 1 ) — > J/ijj + (tttt, rj, n') will be similar in mag- 
nitude, of O(10 — 100) keV. Clearly this is not a rigorous 
result but it does offer a reasonable order of magnitude 
estimate. 

While there are no calculations for radiative transitions 
involving charmonium hybrids there are estimates of ra- 
diative transitions involving hybrids with light quarks 
[37L l38l |. Both calculations found that the El tran- 
sitions between hybrid and conventional states to be 
comparable in magnitude to transitions between con- 
ventional mesons. While neither calculation can be ap- 
plied directly to cc one might take this to suggest that 
the partial widths for ip g (l h ) — > 7 + (J/ip,h c ) and 
ipgiO' 1 , 2 H ) — > 7+ (r/c, Xcj) are the same order of mag- 
nitude as transitions between conventional charmonium 
states. However, it is not at all clear if this extrapola- 
tion to charmonium hybrids is correct as in the flux-tube 
model Close and Dudek [||| showed that the AS = El 
transitions to hybrids only occur for charged particles, 
and hence would vanish for cc. The AS — 1 Ml tran- 



sitions can occur, but are non-leading and less well de- 
fined. Estimates (H for their widths are 0(1 - 100) keV. 
Clearly, given our general lack of understanding of radia- 
tive transitions involving hybrids, the measurement of 
these transitions, ip g — > (00)7, has important implica- 
tions for model builders. 

(Hi) Decays to light hadrons: Decays of the type 
ip g — > light hadrons offer the interesting possibility of 
producing light exotic mesons. Estimates of annihila- 
tion widths to light hadrons will be order of magni- 
tude guesses at best due to uncertainties in wavefunc- 
tion effects and QCD corrections. We estimate the an- 
nihilation widths T[ip g (C = — ) — * light hadrons] and 
r[cc(0 = +) — * light hadrons] by comparing them to 
T(ip' — > light hadrons) and T(i]' c — > light hadrons). 
The light hadron production rate from ip g (C — — ) de- 
cays is suppressed by one power of ct s with respect to 
ip g (C = +) decays. This very naive assumption gives 
?{ip g {C = -) -» light hadrons] - O(100) keV and 
r[cc(C = +) -► light hadrons] - O(10) MeV 
These widths could be smaller because the qq pair in 
hybrids is expected to be separated by a distance of or- 
der 1 / Aqc d resulting in a smaller annihilation rate than 
the S-wave ip' an d Vc states. 



V. SIGNATURES 

The decays discussed above lead to a number of pos- 
sible signals: ip g -> D^>D^**\ ip g (0+-,2 + -) -> J/0 + 
(tt+tt", 77,77'), ip g (l~ + ) -> ?7 C + (tt+tt^, 77, 77'), tpg -> (cc)7, 
and ipg ~ ¥ light hadrons. Of the possible decay modes, 
ip g — > J/i]m + ir~, ipg — > J/ipr), and ip g — > (cc)j give dis- 
tinctive and easily reconstructed signals. In the former 
case the subsequent decay, J/ip — > e + e~ and offers 
a clean tag for the event. 

A good place to search for hybrids in tp g — > J/07r + 7r~ 
is to look for peaks in the invariant mass distributions 
M(e + e~7r~7r + ) — M(e + e~). Babar observed a strong 
signal for the tp' in such a distribution from the decay 
chain B -> ip(2S) + X -> J/-0(e+e~)7r+7r-. Babar's 
efficiency for the 0(e + e~)7r + 7r~ final state is about 20%. 
With 2.3 x 10 7 BB pairs from an integrated luminosity 
of 20.3 ftr 1 Babar observed ~ 972 i}j(2S) -» J/^tt+tt- 
events. Both the Babar and Belle collaborations have 
collected over 100 fb _1 of integrated luminosity and each 
expects to collect 400 fb _1 over the next few years. 

Both the H and 2 H should decay via the ip g — > 
J/ipmr cascade. Although the lattice predictions for the 
H and 2 H masses are above DD** threshold there 
is still considerable uncertainty in these values and the 
flux tube model predicts masses approximately 4.1 GeV. 
We therefore consider both cases, where the ipg lies both 
below and above DD** threshold. For the low mass sce- 
nario, combining our estimates of B(B — > ip g +X) ~ 10~ 3 
and B[ip g (2+-) -> J/^tt+tt-] ~ 0.2 (the 4.1 GeV PSS 
case in Table I) with the PDG value of B(tp -> t + l~) = 
11.81% |40j and the Babar detection efficiency we esti- 
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TABLE I: Decays of exotic Charmonium Hybrids. The fall-apart widths are taken from Page, Swanson and Szczepaniak |32T| 
where IKP refers to the Isgur Kokoski Paton model and PSS to the Page Swanson Szczepaniak model. D**(l~£) and D** 
are the low and high mass 1 + charmed meson states. The high mass state is identified with the Di(2420) state. The hadronic 
cascade decays are rough estimates. See discussion in the text. 



ip g State Final State 


IKP 




PSS 




PSS 








M(V> S ) = 4.1 GeV 


M(tp g ) = 4.1 GeV 


M(t/j g ) = 4.4 GeV 






Width (MeV) BR (% ) Width (MeV) BR 


(%) 


Width (MeV) BR 


(%) 


1"+ 


D*D 


4 


28.4 


0.8 


7.3 


0.1 


0.3 




D**(2 + )D 


- 


- 


- 


- 


0.5 


1.3 




D**{l+)D 


- 


- 


- 


- 


25 


63.5 




D**(l+)D 


- 


- 


- 


- 


3.7 


9.4 






0.1 


0.7 


0.1 


0.9 


0.1 


0.3 




light hadrons 


10 


70.9 


10 


91.7 


10 


39.4 




^Total 


14.1 




10.9 




39 






D*D 


1 


83.3 


0.3 


60 


0.2 


5.1 




D**(2 + )D 










0.5 


12.8 




D**{l+)D 










3 


76.9 




J/ipmv 


0.1 


8.3 


0.1 


20 


0.1 


2.6 




light hadrons 


0.1 


8.3 


0.1 


20 


0.1 


2.6 




^Total 


1.2 




0.5 




3.9 




o+- 


D**(1+)D 










25 


62.2 




D**(1+)D 










15 


37.3 




J/ipmv 


0.1 


50 


0.1 


50 


0.1 


0.2 




light hadrons 


0.1 


50 


0.1 


50 


0.1 


0.2 






0.2 




0.2 




40.2 





mate that for 100 fb 1 of integrated luminosity each ex- 
periment should observe roughly 50 events. If the 2 H 
lies above the DD** threshold the BR for 2+ - -» J/ipmr 
decreases significantly to 2.6 x 10~ 2 (the 4.4 GeV PSS 
case in Table I) lowering the expected number to about 6 
events. Similarly, for the H hybrid we estimate roughly 
1200 events if it lies below threshold but only 5 events 
once the DD** decay modes open up. 

The 1 h state is expected to be the lightest exotic cc 
hybrid [rH l20| and therefore the most likely to lie be- 
low DD** threshold. However, in this case the cascade 
goes to ?7 c 7r7r, a more difficult final state to reconstruct. 
We use our estimate of B(B — > ip g + X) ~ 10~ 3 and 
the estimate given in Table I for the 4.1 GeV PSS case 
of B{ijj g {l-+ -> r/cTT+Tr-) ~ 9 x 10~ 3 . The Babar col- 
laboration studied the decay B — > r/ c K by observing the 
rj c in KKtt and KKKK final states. Combining the 
PDG values for the BR's to these final states with the 
Babar detection efficiencies of roughly 15% and 11% re- 
spectively we estimate that for 100 fb -1 each experiment 
should observe roughly 10 events. If the 1 h lies above 
the DD** threshold, the BR for 1 h — > 7T7r7y c decreases 
to 3 x 10 -3 lowering the expected number to about 3 
events. 

The radiative transition, 4>g(l h ) — > "fJ/ip, also has a 
distinct signal if it has a significant branching ratio. If 
we take the conservative value of T(ip g (l h ) — > jJ/ip) ~ 



1 keV, the BR's for this transition would be rather small. 
On the other hand, a monochramatic photon offers a 
clean tag with a high efficiency. One could look for peaks 
in M(/x + /x~7) — M(/i + /i~). Babar observed Xci and \c2 
this way obtaining ~ 394 Xci's and ~ 1100 Xc2's with 
a 20.3 fb" 1 data sample and an efficiency of about 20 % 
for the J/V>7 final state 0- So although the rate may 
be too small to observe, given the potential payoff, it is 
probably worth the effort to perform this search. Wc 
also note that both Babar and Belle should be able to 
see the 3 Z?2 state via radiative transitions [2£| so that 
even if they do not discover a charmonium hybrid they 
will almost certainly add to our knowledge of quarkonium 
spectroscopy. 

We note that Babar has measured BR's for B — > (cc) + 
h at the 10~ 6 level with 61.6 x 10 6 BB pairs (B(B+ -> 
XcoK+, Xco tt+^-) = (1.46 ±0.35 ±0.12) x Hr 6 ) 
demonstrating the accessability to these levels of com- 
bined BR's. 

Experiments might also look for charmonium hybrids 
in invariant mass distributions of light hadrons. For ex- 
ample, Belle observed the Xco by looking at the invari- 
ant mass distributions from the decays Xco — * 7r + 7r~ and 
Xco — * K + K~ 6]. They found efficiencies of 21% for 
Xco — * 7r + 7r~ and 12.9% for Xco — * K + K~, obtaining 
~ 16 events in the former case and ~ 9 in the latter. 

The decay to charmed mesons also needs to be stud- 
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ied. Because there are more particles in the final state it 
will be more difficult to reconstruct the charmonium hy- 
brid. On the other hand, with sufficient statistics these 
channels will be important for measuring the ip g quan- 
tum numbers and distinguishing their properties from 
conventional cc states. 

The final consideration in charmonium hybrid searches 
is distinguishing the signal from background. The largest 
backgrounds in these final states will be via charmonium 
states with similar masses and fairly narrow total widths 
|28|. The only charmonium states with these proper- 
ties are the missing cc and 1 Dp states whose masses 
are predicted to be ~ 3.8 GeV [42|. They lie below 
DD* and are forbidden to decay into DD because of 
parity conservation. So they are expected to have nar- 
row widths of 300-400 keV and should be easily tagged 
through the dominant El transitions into the Xcj states. 
These states are expected to be produced with branch- 
ing ratios of 0(1%) with B[ip(2—) -> ttttJ/^} ~ 0.12 
^>HEH3' The successes of QCD motivated quark mod- 
els for these states give reasonably reliable predictions 
for their masses so it should be possible to distinguish 
them from the charmonium hybrid on this basis. In any 
case, their discovery would be interesting in their own 
right. For exclusive processes such as B — > ip g + K^* 1 
the If M would have a definite momentum in the B rest 
frame. Careful study of IfM momentum spectra is an- 
other tool that could be used to separate the signal from 
other sources, and seek excess of low momentum K^*' 
recoiling against the ip g ~ 4 — 4.5 GeV. 

VI. CONCLUSIONS 

Establishing the existence of mesons with explicit glu- 
onic degrees of freedom is one of the most important chal- 



lenges in strong interaction physics. As demonstrated by 
the discovery of the rj c {2S) in B decay, B decays offer 
a promising approach to discovering charmonium hybrid 
mesons. In this letter we have described a strategy to 
search for these states. While there is no question that 
our estimates for the various partial widths are crude, 
the essential point is that these states are expected to 
be relatively narrow and that distinctive final states are 
likely to have observable branching ratios. Given how 
much we can learn by finding these states we strongly 
advocate that some effort be devoted to their searches. 



Note added: After the completion of this paper the Belle 
Collaboration published the observation of a new char- 
monium state in exclusive B^ — ► K ± Tr + n~ J/ip decays 
with mass 3871.8 ± 0.7(stat) ± 0.4(syst) MeV which is 
most likely the 3 D C 2 state j43|. However, if this state 
were found to have natural spin parity, + , 1~, 2 + etc, 
then a dynamical suppression, such as expected for hy- 
brids, would be called for. 
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